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Abstract 

In this report, we briefly revisit the formulation of density weighted filtered density function (DW- 
FDF) for large eddy simulation (LES) of turbulent reacting flows, which was proposed by Jaberi et al. 
(Ref. 1). At first, we proceed the traditional derivation of the DW-FDF equations by using the fine 
grained probability density function (FG-PDF), then we explore another way of constructing the DW- 
FDF equations by starting directly from the compressible Navier-Stokes equations. We observe that the 
terms which are unclosed in the traditional DW-FDF equations are now closed in the newly constructed 
DW-FDF equations. This significant difference and its practical impact on the computational simulations 
may deserve further studies. 

1.0 Introduction 

The concept of filtered density function (FDF) method (Ref. 2) was extended to the compressible 
turbulent reacting flow by introducing a filtered mass density function (FMDF) (Ref. 1), which is a 
density weighted FDF and is referred as DW-FDF in this report. The work was considered as the 
extension of the probability density function (PDF) method (Pope (Ref. 3)) to the LES simulation with 
the advantage of its direct simulation of turbulence-chemistry interaction without modeling. 

This report follows our previous work on the conservational PDF equations (Shih and Liu (Ref. 4)) to 
explore different ways of constructing the DW-FDF equations for simulations of filtered turbulent 
reacting flows. It will be shown that there are significant differences between the FMDF equations and 
the newly proposed DW-FDF equations, with the latter requiring much less empirical modeling. 

In Section 2.0, we review the definition of DW-FDF starting from a delta function of a turbulent 
random variable, which is called in the literature as the fine grained probability density function (FG- 
PDF) /'. The definition of a filtered turbulent variable used in LES is also briefly reviewed. Then we 
explore the relationship between the DW-FDF and the filtered turbulent variable. This relationship 
suggests the concept of “DW-FDF mean”, ( ) L which is the analogy to the statistical mean, ( ) in the PDF 
theory. A few properties of the DW-FDF mean are explored, and they are found to be useful in 
constructing the transport equations for the DW-FDF. 

In Section 3.0, we first formulate a compressible transport equation for FG-PDF, and then apply a 
stationary or homogeneous filtering on it to construct the equations that we refer as the traditional DW- 
FDF equations, which are consistent with the FMDF equations introduced by Jaberi et al. (Ref. 1). 

In Section 4.0, we adopt a similar methodology described in Reference 4 to construct the DW-FDF 
equations by directly starting from the filtered compressible Navier-Stokes equations. The resulted 
equations are significantly different from the traditional FMDF equations. For example, the unclosed 
terms appearing on the right hand side of FMDF equations are now closed in the newly constructed DW- 
FDF equations. 
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It is noted here that the quantity DW-FDF (or FMDF) is fundamentally different from the quantity 
PDF, namely, DW-FDF is a random quantity; but PDF is a deterministic quantity. The “DW-FDF mean” 
defines a density-weighted filtered turbulent random variable; but the “PDF mean” defines a turbulent 
mean variable. The DW-FDF, unlike PDF, does not represent a probability density of a random variable; 
it is only a density-weighted filtered FG-PDF and is still a random variable. These fundamental 
differences between DW-FDF and PDF may require a further investigation on their respective solution 
procedures based on the stochastic differential equation methods. 

2.0 Density Weighted Filtered Density Function (DW-FDF) and Filtered 
Turbulent Variables 

In this section, we will first review the basic properties of FG-PDF and its density-weighted filtered 
quantity, DW-FDF, then explore the relationship between DW-FDF and filtered turbulent variable. This 
provides the basis for further establishing the transport equations of DW-FDF. 

2.1 Fine Grained Probability Density Function (FG-PDF) f{j (V ;x,t), /$ (\|/ ; x,t ) 

The FG-PDF for turbulent velocity and scalars (e.g., species mass fraction, internal energy) are 
defined as follows (Pope (Ref. 3)), 

3 

i=i 

N+l 

/®(v; x,O s 8(©(x,o-v) s ]^[5(a ) l -(x,o-Vf) (2) 

/= 1 

where 8 denotes the delta function, U(x,t) is the turbulent velocity vector ( U \ , U 2 , Uj), O(jc,0 is the 
turbulent scalar array (O l5 0 2) - • •, O v , ® w+1 ), for example, N species mass fractions and one internal 
energy cDjv+i = e, the x,t denote the physical space variable (x\, x 2 , x 3 ) and the time t, V = (V\, V 2 , V 3 ) and 
\|/ = ( v|/ , , \|/ 2 , — , \) ) N , '|/ V i i ) are the sample space variables for U(x,t) and 0(.v,t), respectively. The FG-PDF 
has the following properties that will be used in this study: 

i. The Unity Integral of FG-PDF 

I_oo fu{V-,x,t)dV = \, J f^(y\i;x,t)ch\i = ] (3) 

Note that (- oo, + go) represents the whole domain of sample space, for the species mass fractions it really 
means (0,1). 

ii. The Differential Chain 


^ L= ^ L du L= _^ L du L= _^_f atV 

dt dU i dt dVj 8t dV t y U dt y 


(4) 


8xj 8U i 8. Xj 8Vj dxj 


8 


8V ; 


fu 


dUi 

8x; 

J J 


(5) 
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Note that the summation convention on the same index, e.g., “z” is implied. Therefore, we have 




^ + c/ .^0 
v dt ' dx J J 


(for incompressible flow) z, j = 1, 2, 3 


and 


dp/tz dpUjf{j Q 


dt 


dx ; 


dV r - 


fu 


r dpU t dpUjUj ' 


dt 


dxj j 


Note that the continuity equation of compressible flow has been applied here. Similarly, 

ofl _ 50, _ dj% 50, _ 5 ( 50, ' 

7® 


dt 50,- dt 8\\i i dt d\\ij 


5/® _ 5/i 50, _ o/i 50, _ 0 


Gt 


dx; 50, dx : 5v(/, 5x.- ovp , 


50 
/o ~ 

V 8x J J 


and 


5/i +[/ 5/i_ 5 


5t 7 5v, 5\|/, 


/® 


^50, rr 50 ^ 

L + U , 1 

5l 3 dxj , 


V 


(for incompressible flow) ( 


j z = 1, 2, • • • , IV + 1 


U = l,2,3 


5p/i + 5 p ^y7o _ 5 


5l dxj 5\|/, 

iii. The Sifting Property 


/i 


^5pO, 5pO,J7 ^ 


5l 


j_ 

dx j J 


( 6 ) 


(for compressible flow) z, j - 1, 2, 3 (7) 


( 8 ) 


(9) 


( 10 ) 


fi = 1, 2, • • - , TV + 1 , 

, (for compressible flow) ( (11) 

1 7 = 1) 2, 3 


t / (^0/fy( F i jc ,o=t / (^0 6 ( (/ ( jc ’0- F ) =F6 (f / (^0- F )= F /^( F ; jc ,o ( 12 ) 

o (-V, 0/i (V|/ ; AT, o = Y|/ /o (V|/ ; jc, ?) (13) 


More generally, 


Q(u(x,t))f{j{V-x,t) = Q(V)fij(V-x,t) 

Q(®(x,t))f®(\\i ; x,t) = ; *,t) 

where Q( ) denotes a general (but not the time or space differential) function form. 

iv. The Statistical Mean of FG-PDF 

Equations (1) and (2) indicate that FG-PDF is a random quantity, because it is the delta function of 
the random variable U(x,t) or 0(.v,Z). The basic probability theory says that a random variable and its 
(non-differential) function share the same probability density. Therefore, the probability density function 
PDF of the random variable FG-PDF is f v (V; x,t) or x,t), which is the PDF of turbulent velocity 
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U(x,t) or scalars O(x,0, respectively. It is easy to verify that the statistical mean of FG-PDF is the 
corresponding PDF (Pope (Ref. 3)): 


(. fu(V;xj)) = (5(u(x,t)-V)) = \ + _y(V'-V)f u (V';x,t)dV' = f u (V;x,t ) (15) 


Similarly, 


(/<Uv;*>0) = /©(v;*>0 (16) 

The FG-PDF is a useful quantity for mathematical derivations. It is interesting to examine the following 
expressions: 


jV f{j (V ; x,t) dV s {u{x,t)) FG = U(x,t) 

(17) 

J V /© 0 V \x,t) dv/ = (® (x, t)) FG = ® (x, t) 

(18) 


which indicate that the “FG-PDF mean” defines the original random variable itself. More generally, 


jQ(V)f u (V;x,t)dV = (Q(U{x,tj)) FG =Q(U{x,tj) 
J C?(v)/©(V ; x,t)dy = (Q(®(x,t))) FG = Q(<b(x,t)) 


(19) 


2.2 Filtered Turbulent Variables and DW-FDF F v (V ; x,t), /$ (vj/ ; x,t) 

First we review the definition of filtered turbulent variables conventionally used in the turbulent 
simulation and the definition of DW-FDF, and then we explore the relationship between the DW-FDF 
and the filtered turbulent variables. 

2.2.1 Filtered Turbulent Variables 

In the conventional simulation of compressible reacting flows, we often deal with two types of 
filtered turbulent variables: one with the density weighting, the other without the density weighting. The 
filtered turbulent variable without the density weighting is denoted by §(x,t) and is defined as 

— f+CO 

<t>(jc,0= f §(x,t')G(t -t')dt' (20) 

J— 00 

where <)> is the unfiltered turbulent variable, e.g., velocity components U h density p, pressure P, species 

N 

mass fraction ®, and internal energy e = ^ ® m e m . The integration is over the entire time domain 

m = 1 

— oo < t' < + oo. G(t - t") is the time filter with a constant filter width Ay and satisfies the following 
condition and asymptotic property: 


|*+GO 

J— 00 


G(t-t')dt' 


= 1 


( 21 ) 
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p+00 

<t>(jc, t') G(t - 1') dt' = t), 

J— 00 


US A r — ^ 0 


The density-weighted filtered turbulent variable is denoted by (K*M) and defined as 


<i>o,o=^ 


It is noted here that these filtered variables <|)(- V , i ) , (j)(jc,i) represent the large scale turbulent variables 
(Shih and Liu (Refs. 5 and 6)), and they are still random but contain relatively low frequency part of the 
turbulent motion when comparing with the unfiltered turbulent variable (j)(.v,t). 

2.2.2 Definition of DW-FDF 

Jaberi et al. (Ref. 1) introduced the density-weighted space filter operation p(.v\t)G(.v - .*') on the FG- 
PDF f{j(V ; x,t), /([, ( i|i ; x,t) and referred to them as the filtered mass density function FMDF. However, 
we prefer a more intuitive name related to the existing FDF methodology (Colluci et al. (Ref. 2)), and 
refer them as the density weighted FDF (DW-FDF). With the density-weighted time filter p(x,t')G(t - t') 
we define the following marginal DW-FDF as 

Fu(V ; x,t) = J_ + ” p (x,t')f{j{V-x,t') G(t - t') dt' 

= \ +Xj p(x,t')5(u(x,t')-V)G(t-t')dt' 

J “°° (24) 

p+°° . . 

F 0 0 ; x,t) = J ^ p(x,t' )/o(v|/ ; x,t')G{t-t’)dt' 

r +°0 . 

= j p(.\;,?')S(<D(.\;,?')-i|/)G(i -t')dt' 

Obviously, DW-FDF, i.e., F ( { V; x,t) or F'<|,(\|/; x,t), is still a random quantity. It satisfies the following 
“normalization” property: 

00 00 

| Fjjiy ; x,t)dV = | | p(jc,i')/^(K ; x,t')G(t -t')dt' dV = J p(x,t')G(t -t')dt' - p 

-CO -CO ( 2 5 ) 

GO GO 

| (v|/ ; x, t) di// - | | p(jc,i')/^(i|; ; x,t')G(t-t')dt'd\\i-^ p[x,t')G(t -t')dt' = p 


With the definition of DW-FDF described in Equation (24), we can exactly deduce the filtered turbulent 
variables that are conventionally defined by Equations (20) and (23). For example, 

| +C °VF u (V;x,t)dV^j +CO j +OC Vp(x,t')8(u(x,t')-V)G(t-t')dt'dV 

p +00 p+oo . . 

= J J U^xd f )p(xd') 5 ^U^x 9 t')-V)G(t-t , )dt'dV ( 2 ( 

= J p(x,t')U(x,t')G(t -t')dt' = pU = p(x,t)G(x,t) 
where the sifting property of FG-PDF has been used in the second line of (26). Similarly, 
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( 27 ) 


/• +00 ~ 

J V|/F 0 (v|/ ; x,t)d\\i = pO = p(x,?)^>(x,?) 

If we denote the left hand side of Equations (26) and (27) as the operation ( U) L and (0) L , we may write 
Equations (26) and (27) as 


(U) L =j^°VF u (V;x,t)dV = p(x,t)U(x,t) (28) 

(<&) L = J x,t)d\f = p(x,t)<&(x,t) (29) 

and for the function (Q( U(x,t)) or W(0(x,t)), 

(< Q{U)) l = \ + ™Q{V)F v (V ; x,t)dV = p(x,t)Q(U(x,t)) 

(30) 

(W(<S>)) L = j ^(v|/)i+(v|/; x,t)dvf = p(x,t)W(<5>(x,tj) 

Furthermore, we may consider the derivatives VP, VU, VO as the new random quantities and legitimately 
write 


(VP) £ = p VPVP, (VU) L =pVU, (VO) i .=pVO (31) 

However, the operation ( ) L does not have the differential commute property, i.e. 

(VP) L *V(P) L , (VU) L *V(U) L , {V*) L *V{*) L (32) 

because pVP^v(pP), (5 VU ^ v(pf/j, (5VO^V^pO). 

It can be verified that 


I \ C +co _ — ~ — 

{u,Uj) L = ]_ f v t V j F U (V ; x,t)dV = p&tpp j(x,t) (33) 

We can also write ^ Uj&jj for the joint variables as 

(Uj®j ) L = | ^ Vj \\fi F U 0 (1 V , \\i;x,t) dVd\s/ = p(x,t)U J <F l (x, t) (34) 

where, F Ut ®(H,\|/; x,t) is the joint DW-FDF defined as 

Fu t <b (V , \|/ ; x, t) = | p(jt,f)8(i/(jc,f)-F)8(o(jt,f)-i|/)G(f-0<if (35) 

which also satisfies the normalization property: 


00 00 

| j* Fu,<b (V’ V i x > 0 dVd\y = J p(x,t')G(t -t')dt' = p (36) 

—00 —00 
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2.2.3 DW-FDF Mean < ) L 

Analogous to JVfj (V; x,t)dV = (U which defines the statistical mean of the random velocity 
by its PDF f v , we may refer the operation < ) L as the “DW-FDF-mean”, which defines the density 

weighted filtered turbulent variable, e.g., (u) L = p U , (Uj®^ = pt/ ; 0 ; - (see Eqs. (28) to (34)). From 
Equation (35), we may follow Pope (Ref. 3) to define a “conditional” DW-FDF on the condition O = \|/ as 


F U\^ V \^'’ X F) 


F U(S? (V,\\i-x,t) 

F<b(y;x,t) 


(37) 


and the “conditional DW-FDF-mean” as 


(U (x, t)\ \\i) L = VF u \ <b (V\y,x,t)dV = — - J* V F U ([) (V, y;x,t)dV 

1 f+OO /•+00 . . . . 

= — j J Vp(x,t')h(u(x,t')-V)h(Q>(x,t')-y)G(t-t')dt'dV (38) 

At J_co _GO 

= — | p(.v,f)t/(jt,f)S(o(x,f )- \\i)G(t -t')dt' 

Fq j - co 


Then we have 


j ^ P (*, t') U (x, f )/o ('I' ; ?)G(t - 1') dt' = | ^ V F Uf(b (V, i|/ ; x, t) dV = F 0 ■ (u (x, t) | v| i) L (39) 


and 


(•+00/ , i \ \ C+GO C+GO 

J_oo ( F $ ' i U ( x ’ v ) L ) d V = }_„ V F U , ® (V,y,x,t) dVdvf 

C+00 . 

= j p(x,t')U(x,t')G(t-t')dt' 

-p(x,t)U(x,t) 

Equation (39) can be extended to any other turbulent quantities, for example, VP, Vt/, VO, 5,(0): 


(40) 


r +°° / \ 

\_ 7 P (x,f ) VP(x, f) f{j ( V ; *, f) G(t - f) dt' = F v • (VP| v) L 

J_ + J p (x,t')VU(x, t') f{j (V; x, t') G(t - t') dt' = F v • (’ VU\ v) L 

] p(x,t')VQ>(x,t') x,t')G{t -t')dt' = F® ■ L 

P (x, f ) 5, (0(x, o) /o ( V; A t') G(t - f) dt' = P 0 • (Sj (0)1 v) L =F <b -S i (y) 
where VP, VU, VO, are viewed as the random variables in addition to P, U, O. 


(41) 
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3.0 Traditional DW-FDF Equations for Fjj ( V ; x,t),F ( ^ (\| / ; x,t) 

In this section, we briefly review the traditional way of deriving DW-FDF equations starting from the 
equation of FG-PDF. The “exact” but not closed DW-FDF equations are derived, which are consistent 
with the traditional FMDF equations introduced by Jaberi et al. (Ref. 1). 

3.1 Traditional DW-FDF Equation for F v (V ; x,t) 

The transport equation for the fine grained PDF with the variable density p(x,f) can be written as (see 
Equation (7)): 


dpfu [ ^ U jfu 
dt dxj 



f 


f 


_dP + d 

v dXf dxj v 


2pS, 


2 A 
~^\ iF> ij S kk 


(42) 


Where the right hand side of (42) is from the Navier-Stoke equation (47). 

Applying the time filtering operation G(t - t") with a constant filter width on Equation (42) and using 
the sifting property of FG-PDF, we may obtain (noting that the differential commutation property of the 
time filtering operation is valid with a constant time filter width) 


dF u , v dF u _ d 

dt i dx: dV: 

J 1 l -CO 


f p f{j ~ G(t-t')dt' + f p f{j 

J n nr. J 


vP dx u 


1 8 


v p dXj 


2p5,--pS,5 ( 


3 


ij°kk 


\G(t-t')dt' 


According to Equation (41), the DW-FDF equation for the velocity can be written as 


8F, 


u 


dt 


+ V; 


dF, 


u 


J 8x j 


dV: 


-F, 


ir\ 


1 8P 

p dxj 


V) +F 


1 8 


u-\ 


p dx 


7'V 


-—^ijSkk 
4 J 


V 


(43) 


At this point, the velocity DW-FDF equation is general but unclosed because of the unknown terms of the 
conditional DW-FDF means. The last term in Equation (43) is referred as the molecular mixing term and 
was modeled in various empirical ways by different researchers. Later, we will see that, when using a 
different way to derive the DW-FDF equation, the molecular mixing term is automatically closed with no 
need of modeling. The same happens to the pressure gradient term. 


3.2 Traditional DW-FDF Equation foriy t) (\|/ ; x ,t) 

Applying the time filtering G(t - t') on the following FG-PDF equation for scalars, 


dp/<p [ d 9 u jf<b 5_, ,, 

8t dxj S\\ij ° 


f d f 


\ dx j 


P r (0 


dx i J 


+ p5,-(®(x,0) 


(44) 


we obtain the DW-FDF equation for the turbulent species: 


8F, 


$ 


dt 


+ - 


d [ F a>( u j V 


dXj 




r, I 1 5 50, 

F ® '( _ 7^P r (0 


p dx 


J 


dx 


v|/) +F <b -S i (y) 

L 


(45) 
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where -SVn(v|/) = 0- Equation (45) is also unclosed, because the conditional DW-FDF means need to be 
modeled. Later, we will also see that using a different derivation of DW-FDF equation, the molecular 
mixing term is closed. 


4.0 DW-FDF Equations Derived From Filtered Navier-Stokes Equations 

4.1 Filtered Compressible Navier-Stokes Equations 

The Navier-Stokes equations for a compressible reacting flow can be written as 


gp | 8 P U J _Q 

dt dxj 


( 46 ) 


op Uj d p UjU j _ dP d 
8t dxj dxj dxj V 

dpe dpU i e_ dq t 


2 P 


f 


dt dx; 


dx. 


+ PS k k +2pv 


1 




SjSj—SuS^ 
v 3 


+ Q 


j 


5pO„, dpU^ h 


dt 


- + 


dxj dxj y 


pr 


( m ) 

dX: 


+ w m m = \,2,-,N 


N 


p = pR^ t 


n = 1 W n 


-pKC^-fp 

OX: . OX: 

1 m = 1 1 


(47) 

(48) 

(49) 

(50) 

(51) 


Applying the time filtering operation with a constant filter width on the above equations, we obtain the 
following filtered Navier-Stokes equations: 


dp | dp Uj _ q 
dt dxj 


dpUj dp UjU j 
dt dxj 


_dP d 

dxj dxj 


2 P v (%-^8y%) 
V 4 


dpe dpUje 
dt dxf 


-~^- + PS kk + 2pv 

ox. 


Si i S ii --S ii S kk 


+ Q 


gpjm , dp Uj® m 

dt dxj 


dx , 


pr 


(/») d® m 
dx , 


\ + W n 


i J 


m = \,2,---,N 


(52) 

(53) 

(54) 

(55) 


where 
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1 w„ c., , w„ 

n = 1 ' n = 1 « 


or 


p = y 1 ^^tv+i 

c v n = 1 


_ 37 ^ p(/»)/ 

?,- = -pKc 0 - — h » 


dx, 


1 m= 1 


&,• 


( 56 ) 


( 57 ) 


In the above equations, k, v and T ( "° are the molecular heat conductivity, kinematic viscosity and the 772-th 
species diffusivity, they have the same dimension (i.e., velocity • length). The h m , T are the enthalpy of 
species and temperature, Q is the radiation rate, W m = p.S'„, is the chemical generation rate of the 772-th 
species, Ojv+i represents the internal energy e, R is the universal gas constant. These equations are general; 
however, unlike the constant density flows, further approximations for the terms on their right hand side 
are required in order to complete the filtering process. One of such approximations leads to 


2 pv 


S, , — 8 ,,S, 


ij^kk 


dUi dUj 2 8 U k 
— - + — - — 8 ,,' — - 
8 . %j dxj 3 7 8 x k 


( 58 ) 


In which, we have basically neglected the variations of p and pduring the filtering process, the value of p 
will be considered as the function of P, f, ■ • ■ . Another type of approximation writes 


Similarly, 


2 pv 


S;; — 8 ,,S. 


ij°kk 


'dpUj , dpUj 2 g dp U k ' 
v dxj dxj 3 9 dx k ^ 


<li=-pKc v 


8T_ 

8X; 


N 


Xpr < ”’\ 


m = 1 


8 X: 


! - KC, 


dpf 

dx t 


N 


Zr 

m = 1 


O) 


dx f 


( 59 ) 


( 60 ) 


pr 


0 m ) 


dx; 


-.(m) ®rn 
8 x ; 


Furthermore, invoking the turbulent kinetic energy dissipation rate: 


2 pv 


SySy-^uS^ 


= ps 


Where v, k, c v and T 1 " 0 are considered as the function of P,f , • • • . 

Therefore, the filtered Navier-Stokes equations can approximately be written as 


dp | dp Uj _ q 
dt dxj 


( 61 ) 


( 62 ) 


( 63 ) 
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(64) 


dp Uj | dpU t Uj _ 8P | d 


dt 
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V 

A 


Ox,- 
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i J 
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or 


A 


- Ql< v- '[>„'[> 




w^A+l 


C v ^ W n 


Vi = - 


dp f 

dx, 


A 


1 m = 1 


y r (m) dp h m Q>m 
^ dx ; - 


(65) 

( 66 ) 


(67) 


( 68 ) 


These equations are still considered as quite general, because i) they are exact if the flow becomes 
incompressible, ii) all the approximations made in Equations (64), (65) and (66) are related only to the 
molecular diffusion terms that are less important and even negligible comparing with the convection 
terms on the left hand side for turbulent flows at high Reynolds number (see Tennekes and Lumley 
(Ref. 7) and Pope (Ref. 3)). For the LES simulation, Equation (63) to (66) are often used together with the 
further approximations for (67) and (68): 


A 




1 W ,1 

n = 1 n 


f 

V 

v AT 


pRT^^RT J g = y® iL 

\/f * J in 


A 


Vi = - KC\ 


M 

dpf 

8X; 


n = 1 W n 


(69) 

(70) 


The momentum flux p £/,-£/,- , the energy flux p Up and the species flux p(/,® ;;; are considered to be 

critically important in LES simulations and should be carefully modeled. Many models in the literature, 
from the simplest Smagorinsky (Ref. 8) model to the complex two-equation models (Refs. 9 and 6) 
including the dynamic procedure (Ref. 10) have been suggested. In the next Section 4.2 and 4.3, we will 
derive the DW-FDF equations directly from Equations (63) to (68). 


4.2 DW-FDF Equation for Fp(F;x,f) 

Using Equations (28) and (33), the left hand side of Equation (64) can be written as 


dpUj dpUjUj 
dt dxj 


h\ v ‘ F " dv+ l:\ v ‘ v i F " 

—00 J —00 



—00 V 


dFy ' 

dx j 
J 


dV 


(71) 
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The pressure gradient term can be written as using (67) 


dP 

dxj 


dx ; 


^ P R 

V C v n = 1 


n^N + 1 



( 72 ) 


During the above arrangement, no approximations have been made other than the integration by parts and 
a zero integration property like Eq.(74). 

Note that the pressure gradient term can be closed only when the joint DW-FDF, F UA> , is considered, 
otherwise, the “conditional DW-FDF mean” will be unavoidable. Now the molecular mixing term in 
Equation (64) can be written as 


dx ; 


dp U t dpUj 2 dp U, 


— 8 „- 


. dxj dxj 3 lJ dx m 


dx i 


\ 8x j 


V t F v dV + 


dx,- 


2 d 

V j F u dV — 5,. 

7 u 3 J dx„ 
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V m F v dV 


V 
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—00 
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dx k 

V K -00 
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3 J dx, 
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dxj 


d_ 

dx 


V: 


—00 

or 

00 a 

-k— 

J dx.- 

—00 J 


V J —00 

f d dF v ' 

y dXj dxj j 


: J j v^v^yv- h tJ1 F ] v,F-(v m F u pv 


d_ 

dx k 


dVr 


3 J dx n 


_d_ 
dVk 


(73) 
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l v - 


d 

dxj 


f d 2 Vj Fjj ^ 
v 


V 


dV k dx k 


I 


( a2 


dV+ I Vi-hj — 
3 dx J 


dXnFu 

dV k dx m 


dV 


( 


dVkFu 

dV k dx 


f ^2 


J J 


dV - V: 

J dx,- 

—00 J 


d Vj Fjj 

dV k dx k 


( *2 


dV+ J 


dV m Fu 

dV k dx m 


dV 


Note, on the third line of (73), we have applied the integration by parts and the following type of zero 
integration similar to the one used in PDF formulations (Ref. 3): 


J —{V i V k F u )dV = 0 (74) 

—CO k 

Finally, collecting the terms that are in the integrands and factored by V h we obtain the transport equation 
for F,, ( V; x,t) as 
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(75) 


In Equation (75), all the terms are closed, except for the one involving the conditional DW-FDF mean 
originating from the pressure gradient. It should be noticed that one of the most important terms in the 

momentum equation, i.e., p U-Uj , now is closed in the corresponding DW-FDF equation. And the less 
important molecular diffusion terms remain in the closed form. 


4.3 DW-FDF Equation for ;x,t) 

Similarly, we may obtain the DW-FDF equations for the scalars (i.e. species and internal energy) 
from Equations (65), (66) as follows: first, we write the terms on the left hand side of Equation (66) as 


8 f f 8F® 

. ~ pr Vm F <b d V = I 

dt dt J J dt 


(76) 
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dxj dxj 
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(77) 


Then, the terms on the right hand side of (66) can be written as 
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(78) 
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pS m = { S m (w)F 0 dy = - J 


9{S k F c ,) j 

dVk 


(79) 


Where in Equations (78) and (79), we have applied the integration by parts and zero integrations similar 
to (74). Collecting all the integrand terms that factored by \| r m , we obtain 


dt 

or 



A)_ 1 
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V(J n) dF © ' 
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8x i 
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l dx i J 
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(80) 
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dt 
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+ F 0 - 5* (\|/ )L m = l,2,---,N + l 


i j 


This equation also represents the equation of internal energy m = N + 1, where S N+l (\\/) = 0 and other 
source terms in Equation (65) are neglected as suggested by Jaberi et al. (Ref. 1). Equation (80) is 
essentially closed if we consider the joint DW-FDF F UA ,. For the marginal i>, the convection term is not 
closed because of the conditional DW-FDF mean • Then, this critically important term, 

corresponding to pt/ ( 0 ? in Equation (66), must be carefully modeled while the less important molecular 
diffusion term remains in the closed form. In addition, we noticed that the equally important chemistry 
source term pS^lO) in Equation (66) contains complex processes of turbulence-chemistry interaction, 
which is considered as very hard to be modeled. However, it is closed in the DW-FDF equation with no 
need of modeling. This direct calculation of turbulence-chemistry interaction represents one of the unique 
features of DW-FDF methodology. 


4.4 Modeling of Unclosed Terms 

Consider the approximation described by Equation (69), i.e. P ~pRT / M , we may write from 
Equation (72) 


8P __8 
dxj dxj 
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dV k dx k 

The last step in Equation (81) employs a rough approximation: 


( 81 ) 
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(82) 


(°jv+i| f ) z *®n+i =e 

As for the convection term in the species Equation (80), , we may start from a more general 

model for the term p[/ ; O m (see Ref. 6): 


pUp m =P U^ m -T ( T m) ^L 
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- r^” * ~^[c\ Sy + c 2 ciij ) 


opO, 
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(83) 


Where c\ = C2 = - 0.24. This will lead to the following model by directly applying the Equations (27) and 
(39) described in Sections 2.2.2 and 2.2.3: 
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4.5 Summary 

With the models given by Equations (8 1 ) and (84), the marginal D W -FDF equations for F v { V; x,t) 
and E®(v|/; x,t) can be written as 
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It can be verified that the DW-FDF Equations (85) and (86) can exactly deduce the filtered compressible 
Navier-Stoke equations (64) to (66) with the assumed approximations of (81) and (83). However, by no 
means, the models described by Equations (81) and (84) are unique. Furthermore, the variables ( O v+] , 

Uj , Sjj and Q ; y ) are considered as available during the solution procedure of the DW-FDF equations. 


5.0 Concluding Remarks 

We have revisited the derivations of the traditional DW-FDF or FMDF equations by starting from the 
equation of FG-PDF that contains all the terms on the right hand side of the Navier-Stokes equations, see 
Equations (42) and (44). The resulting FMDF equations contain the “conditional DW-FDF mean” 
quantities for all the terms on the right hand side (except for the reaction term), see Equations (43), (45), 
which need empirical models to make the equations closed. 

From the relationship between the DW-FDF and the filtered turbulent variables, it is possible to 
construct the DW-FDF equations directly from the filtered compressible Navier-Stokes equations. Such 
DW-FDF equations have an outstanding feature that all the terms that were not closed on the right hand 
side of traditional FMDF equations are now in the closed form, see Equations (75) and (80), except for 
the pressure gradient term that involves joint DW-FDF. 

In the marginal DW-FDF equations, the two unclosed terms, one related to the pressure gradient and 
the other related to the species convection, are expected to be more important than the molecular “mixing” 
or “diffusion” terms, especially in the case of high Reynolds number. Their physics-based models, such as 
Equations (81) and (84), should be further developed and evaluated. 

The fundamental differences between the DW-FDF (a random quantity) and the PDF (a deterministic 
quantity) warrant further investigation on their respective solution procedure using the stochastic 
differential equation methods. 
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